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Abstract 

We examine the scaling of geodesic correlation functions in two-dimensional 
gravity and in spin systems coupled to gravity. The numerical data support 
the scaling hypothesis and indicate that the quantum geometry develops a non- 
perturbative length scale. The existence of this length scale allows us to extract 
a Hausdorff dimension. In the case of pure gravity we find dn ~ 3.8, in support 
of recent theoretical calculations that dn — 4. We also discuss the back-reaction 
of matter on the geometry. 



1 Introduction 



Remarkable strides have been made in recent years in our understanding of the 
properties of two-dimensional quantum gravity |]^. Calculations carried out within 
the framework of conformal field theory have yielded the gravitational dressing of 
integrated matter field operators, correlation functions on the sphere and the torus 
partition function. On the other hand matrix models have provided us with a pow- 
erful calculational tool that enables us to compute the above mentioned quantities 
and also to perform the non-perturbative sum over topologies. 

Nevertheless there are still important geometrical quantities of physical in- 
terest that are not well understood analytically. Perhaps the most fundamental is 
the intrinsic Hausdorff dimension of the typical surface generated by the coupling of 
matter to 2(i-gravity ||2|, ^ . One may think of the Hausdorff dimension as an order 
parameter characterizing possible phases of the theory. If there exists a power law 
relation between two reparametrization invariant quantities with the dimension of 
volume (y) and length (L), this provides a well-defined fractal dimension (dn) via 
V oc L'^" . As there is no natural notion of a length scale in these theories, one has 
to be introduced by hand, at least in the continuum formulation. In the discretized 
approach this length scale is provided by the short distance cut-off corresponding 
to the finite elementary link length. 

Recently a transfer matrix formalism utilizing matrix model amplitudes has 
been developed that predicts the Hausdorff dimension dn = 4 for pure 2d gravity 
1^]. This approach has not yet been extended to the case of unitary minimal models 
coupled to gravity. On the other hand the analysis of the diffusion equation for a 
random walk on the ensemble of 2d manifolds determined by the Liouville action 
yields a prediction for the Hausdorff dimension which agrees with the transfer matrix 
approach for pure gravity. It may also be extended to include the coupling of 
conformal matter of central charge c < 1 [^. 

These analytic predictions for the Hausdorff dimension rely on the validity 
of certain scaling assumptions. It also appears that there are several potentially 
inequivalent definitions of an appropriate fractal dimensionality. It seems very 
worthwhile therefore to explore these issues numerically. Earlier numerical work 
addressing this question has been remarkably inconclusive 0, Indeed for a 
while it was claimed that there was no well-defined Hausdorff dimension in the case 
of pure gravity In contrast clear numerical evidence for a fractal scaling of 
gravity coupled to c = — 2 matter was found in ||9|. 

In this letter we establish numerically that this scaling behavior is valid for 
pure gravity as well as the Ising and 3-state Potts models coupled to gravity. We 
employ a careful finite size scaling analysis of appropriate correlation functions. For 
pure gravity we find dn ~ 3.8 in qualitative agreement with [^, |5|. For the Ising 
and 3-state Potts models the values of dn that we obtain do not seem to detect the 
back reaction of matter on the geometry. 
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This paper is organized as follows. In section 2 we describe the application 
of finite size scaling to loop-loop correlation functions. In section 3 we outline our 
numerical procedures and results. In section 4 we present the existing theoretical 
predictions for the Hausdorff dimension. Finally section 5 is a discussion of our 
conclusions. 

2 Scaling 

Finite size scaling is a well-established technique for analyzing the critical behavior 
of conventional statistical mechanical models In numerical studies of quantum 
gravity it has traditionally been employed in a rather limited context - typically by 
extracting a power law scaling for integrated matter field operators at the critical 
point. 

In general, the scaling hypothesis asserts that near a critical point an ob- 
servable O, a function of two variables x and y, will depend on only one scaling 
combination /x = y/x^ up to an overall power factor x^ 

0{x,y)r.xPf{y/x'^). (1) 

The powers p and q are related to the critical exponents of the model. We will test 
this hypothesis by analyzing geodesic correlators defined on dynamical triangula- 
tions. 

The fundamental objects in two-dimensional gravity are loop-loop correla- 
tors. To define these consider two marked loops of length I and /' on a triangulation. 
If we define a geodesic distance r between the loops on the graph as the minimal 
number of links that must be traversed to go from I to I', we can define a correlation 
function ni^i/ (r) as 

ni^vir,N)= ^1. (2) 

TeT2{N) 

In this expression T2{N) refers to the class of triangulations with triangles and 
two loops of length I and I' separated by a geodesic distance r. As defined above 
ni^i'{r,N) is proportional to the number of triangulations satisfying the above con- 
straints. We chose to work in the micro canonical ensemble as it is convenient com- 
putationally and the effect of restricting to fixed volume can be exploited in the 
finite size scaling analysis. The point-point correlator n(r,N), which counts the 
number of triangulations with two marked points separated by geodesic distance r, 
can now be obtained from Eq. (|2|) in the limit that the lengths of the loops are taken 
to zero. 

The scaling hypothesis applied to n (r, A^) implies 

n{r,N) = NP n{r/N'i) . (3) 

The combination Iq = N'^ constitutes a dynamical length scale which appears non- 
perturbatively in the theory. It can be used to define a Hausdorff or fractal dimension 
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dn = 1/9 characterizing the quantum geometry. Notice that in this case the expo- 
nent p is not free — it is constrained by the fact that the integral of n (r, A^) over all 
geodesic distances recovers the total number of points N. This yields p = l — l/dn- 
It is easy to measure the point-point correlator numerically and thus determine du- 
This discussion can be generalized to include spin models coupled to gravity. 
In this case the boundary loops will be dressed with fixed boundary spin configura- 
tions. For the point-point correlator we can use the symmetry of the spin models 
to reduce the possible correlators to two distinct types, which we denote /i (r, A'') 
and /2 (r, N). The correlator fi counts the number of points at distance r for which 
the spin variable is in the same state as the initial marked point. The correlator /2 
counts the number of spins in a different state from the initial marked point. The 
total number of points at geodesic distance r is then 

n{r,N)=h{r,N)+f2{r,N). (4) 

The scaling hypothesis can be applied as before to these correlators, resulting in a 
definition of dn- 

To define spin-spin correlators we note that the spin variables of the g-state 
Potts model may be taken to be the unit vectors Cq, of a hyper-tetrahedron in g — 1 
dimensional space. We can then identify the product of two spins as the scalar 
product of the associated link vectors 

aa(Ti3 = Ca- e/3 = ^1 + ^"^y^ <^a/3 " '~~[- 

The (unnormalized) spin-spin correlator with one marked point i is 

gunir, iV) = E E E ^^^J ^id,j - r) e'^^i''^ ^) . (6) 

T {<7,} j 

In terms of the distributions /i and /2 defined earlier we see from Eq. (|5|) that this 
may be written as 

9un{r, N) = h {r, N) - (r, N) . (7) 

q-1 

The scaling hypothesis for gun takes the form 

gun{r,N) = N^''gun{r/N'). (8) 

The overall power is again determined from the constraint that the integral of 
guni'T^N) is just the usual magnetic susceptibility, which scales at criticality as 

X ~ N'"^H , If we make the assumption that these critical systems contain only one 
length scale then it is natural to assume that both n{r, N) and gunir, N) depend on 
the same scaling variable. This implies that s = l/dn, where the Hausdorff dimen- 
sion dn is now that appropriate to the matter-coupled theory. We shall re-examine 
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this assumption critically in light of the numerical results in the final section of the 
paper. 

We will also consider the normalized spin-spin correlator with one marked 
point i: 

9n{r, iV) = E V - ^ ^ - • 

3 Numerical Simulations 

To investigate the validity of the scaling hypothesis we have performed Monte Carlo 
simulations on three models; pure gravity (central charge c = 0), the Ising model 
(c = 1/2) and 3-state Potts model (c = 4/5) coupled to gravity. In the microcanon- 
ical ensemble the partition function of these models is given by 

z{p,n)=J2Zm{P,n) (10) 

TeT 

where Zm{N,(3) describes the matter sector (absent for pure gravity). For a g-state 
Potts model this is 

Zm(/3, iV) = E exp I /3 E (5,,,,^. - 1) ) , (11) 

{or} \ <i,j> J 

were crj denote the Potts spins, i denotes a lattice site and < i,j > indicates that 
the sum is over nearest- neighbor pairs on the lattice. 

The integration over manifolds is implemented as a sum over an appropri- 
ate class of triangulations T. Since it has been observed that finite size effects in 
numerical determinations of critical exponents are generally smaller if one includes 
degenerate triangulations in T, i.e. triangulations allowing two vertices connected 
by more than one link and vertices connected to itself [jll|Q we will work in that 
ensemble. 

In the simulations a standard link-flip algorithm was used to explore the 
space of triangulations and a Swendsen-Wang cluster algorithm employed for the 
spin updates. Lattice sizes ranging from 500 to 32000 triangles were studied and 
typically 10^ to 4 x 10*^ Monte Carlo sweeps performed for each lattice size (a sweep 
consists of flipping about N links and one SW update of the spin configuration). 

3.1 Pure Gravity 

We start with the results for pure gravity. Here we measured the point-point distri- 
butions n(r, A^) both on the direct and the dual lattice. On the dual lattice geodesic 

^This corresponds to allowing tadpoles and self-energy diagrams in the dual lattice formulation. 
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Figure 1: Volume scaling of (a) the location of the peak tq in the distributions 
n(r, N) and (6) their maximal value n(ro) in the case of pure gravity. Data is shown 
both for the direct and dual lattices and the extracted values of dn are included. 

distances are measured as shortest paths going from one vertex to another. Having 
measured these distributions for different lattice sizes there are several ways we can 
use the scaling assumption Eq. @ to extract dn. We use two methods. 

First we fitted a distribution (for a given lattice size) to an appropriately 
chosen function from which we located the maximum of the distribution tq and its 
maximal value n(ro). Then the scaling assumption implies that ro ~ N^/'^h a,nd 
n(ro) ~ N^-^/'^H . We fit to the function 

P/(r) exp(-ar*) . (12) 

The exponential is included in order to capture the long-distance behavior of the 
distribution and Pi is an /-order polynomial. The order of the polynomial is chosen 
in such way that we get a reasonably good fit; a Ath order polynomial turned out to 
be sufficient. We checked that the values of ro and n(ro) did not change appreciably 
if we increased the order of Pi{r). The values of tq and n(ro) obtained in this way 
are plotted in Figs, la and 1& on log-log plots. As expected both quantities scale 
well with (significantly better for the direct lattice). The Hausdorff dimensions 
extracted from the slopes are listed in Table 1. 

Another way to extract dn is to use the scaling relations directly to collapse 
distributions for different lattices sizes on to the same curve, using only a single 
scaling parameter dn- This we have done including all the data (for A^ > 1000) 
and also, to explore the finite size corrections, using only pairs of datasets (A^ and 
2A^). The same functional form Eq. ([T2| ) was used in the fits. The results are shown 
in Table 1, together with the quality of the fits (x^ = / dof). The errors quoted 
indicate where changes by one unit from its minimal value. In Figs. 2a and 2h 
we show overall scaling plots for n(r, A^), both for the direct and dual lattices. 

From these results we can immediately draw a number of conclusions. Con- 
sider first the direct lattice. Fig. 2a shows that the scaling hypothesis is indeed well 
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Direct lattice 
dH 5d 


Dual lattice 


(a) 

126 - 250 
250 - 500 
500 - 1000 
1000 - 2000 
2000 - 4000 
4000 - 8000 
8000 - 16000 
16000 - 32000 


3.640(60) 44.6 
3.707(45) 13.0 
3.727(42) 8.0 
3.770(38) 4.2 
3.800(54) 2.3 
3.804(55) 1.5 
3.810(55) 0.97 
3.830(50) 1.4 


2.497(37) 49.2 
2.715(40) 29.1 
2.871(38) 20.5 
2.996(26) 22.6 
3.111(39) 12.5 
3.217(47) 9.7 
3.264(34) 6.9 
3.411(89) 4.8 


(h) 

1000 - 32000 


3.790(30) 13.0 


3.150(31) 85 


(c) 

position 
height 


3.835(59) 0.03 
4.040(98) 0.09 


3.133(43) 10.45 
3.594(77) 0.37 



Table 1: Extracted values of dn from n(r, in the case of pure gravity. The values 
in (a) are obtained by collapsing data for two consecutive lattices sizes on a single 
curve using one scaling parameter, (h) is the same except data from all lattice sizes 
between 1000 and 32000 triangles are used. In ('cj the values are obtained from the 
volume scaling of ro and n(ro) separately. The quality of the fit is indicated by 
and the errors m. (a) and (h) indicate where changes by one unit from its minimal 
value. 




Figure 2: Scaling plots for the point-point distributions n{r,N) in the case of 
pure gravity; (a) the direct and (b) dual lattice. Shown are the curves fitted to 
distributions after rescaling. The value of dfj is chosen so as it minimized the total 
chi-square of the fits. 
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satisfied for tlie distribution n(r, A^). This is also evident from the low values of \^ 
for the fits (Table 1). The values of du obtained from the scaling of ro and n(ro) 
and also from collapsing the data are close to the expected value of du = 4. These 
results are obtained on moderately small lattices, illustrating the superiority of this 
method of extracting du to earlier numerical attempts. 

But we also notice that there is a systematic increase in the value of du with 
lattice size. Even though this effect is too small compared to the uncertainty in the 
measured values to allow reliable extrapolation to infinite volume du^ it indicates 
that the difference between measured and expected values of du is due to finite-size 
effects. The improvement of the values of the fits with increasing lattice size also 
implies diminishing deviations from scaling. 

It is also intriguing that the scaling of the peak heights seems to give better 
values of du (close to the theoretical results for the direct lattice) . Since the heights 
of the peaks take continuous values, as opposed to the discrete geodesic distance, it 
plausible that they are less sensitive to the discretization 

On the dual lattice we observe much larger finite size deviations. This is 
evident both from Fig. 26 and the values of in Table 1. This is not hard to 
understand. The short distance behavior of n(r, A'^) is dominated by a power growth 
r'^H-'^ _ But as the order of vertices on the dual lattice is fixed to be three, the growth 
of n(r, is bounded by the function 3 x . If du = 4 this means that for small 
values of r the distribution n(r, A^) may not grow fast enough to display the correct 
fractal structure. Only when the lattices are big enough so that the first few steps 
are negligible can the dual lattice be used to extract du- This constraint on the 
growth is not present on the direct lattice, which is thus better suited for extracting 
du- 



3.2 Coupling to matter 

To see how the point-point distributions (and du) change as we include coupling 
to matter we looked at both the Ising and 3-state Potts models coupled to gravity. 
These models are chosen because in both cases the exact solution of the models is 
known^ knowing the exact critical coupling makes the simulations much easier. 

As shown in the case of pure gravity it is preferable to measure on the direct 
lattice and so we have placed the spins on the vertices. In that case the critical 

■^The 3-state Potts model coupled to gravity has just recently been solved using matrix model 
techniques The numerical simulations we do here verify that the solution is correct. To obtain 
the critical coupling from one has to do some reformulation. This leads to /3* — 1/2 log [(45 — 
•\/(45))/(-^/(47) — 2)]. This is for the spins placed on triangles. To get the coupling for spins on 
vertices we use the duality transformation for the g-state Potts model (e^''^ — l) {(P'^" — l) = g [Q. 
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Ising model 


3-state Potts model 


Exponent 


Measured 


Exact 


Measured 


Exact 




0.167(3) 


1/6 


0.199(4) 


1/5 




0.653(8) 


2/3 


0.608(6) 


3/5 




0.318(12) 


1/3 


0.382(30) 


2/5 



Table 2: Comparing critical exponents, obtained using finite sizes scaling in /3c, to 
exact values, for the Ising and 3-state Potts models coupled to gravity 




1000 10000 1000 10000 

iV N 



Figure 3: Volume scaling for ro and n(ro) for the distributions wc measured for 
Ising model coupled to gravity. The same scaling behavior is used to extract dn 
from the slope as in the case of pure gravity, except for gun{^o)- There we used 
n(ro) ~ ^/fdn — l/c^if, substituting the exact values for ^/udn- 



couplings are (as we include degenerate triangulations) : 



13c 



13 + V?' 
14- V7 



(Ising) and (5c 



41 + v/47' 
47-2^47 



(3 — state Potts). 



(13) 

To verify that these are indeed the correct couplings we have performed a stan- 
dard finite size scaling analysis of some observables related to the spin models; the 
average magnetization M. ~ N~^/'^'^" , the magnetic susceptibility % ~ j^'i/^^Ih ^ 
and the derivative of Binders cumulant dBC /dj3 ~ N^/'^'^h _ 'pj^g measured critical 
exponents are shown in Table 2, together with the exact values with which they 
agree very well. The main reason is, of course, that we know /3c, but also including 
degenerate triangulations and placing the spins on vertices reduces finite-size effects 
dramatically. 

Now to the distribution functions. The placement of spins on the vertices al- 
lows us to measure several combinations of distributions; /i(r, N), /2(r, N), n{r, N) 
and gun{f,N)- We have analyzed these distributions in the same way as for pure 
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Figure 4: Collapsing the data for n{r, N) and gun{f^ on a single curve using one 
scaling parameter in the case of an Ising model coupled to gravity. 



gravity. In Figs. 3a and h we show the scaling with volume of rg and n(ro), obtained 
from fitting the distributions to the functional form Eq. (|l2|) . These plots are for the 
Ising model but plots for the 3-state Potts model are very similar. The extracted 
Hausdorff dimensions, for n(r, N) and gun{T, N), are shown in Table 3. As for pure 
gravity we also scaled all the data (for N > 1000), and for pairs of distributions, on 
a single curve. Resulting optimal values of dn are listed in Table 3. The quality of 
the scaling is shown in Figs. 4a and b, where we show scaling plots for n(r, N) and 
gun{r, N) (for the Ising model). Again the value of dn that minimizes is used to 
scale the data. 

In the case of the spin models we also measured the normalized spin-spin 
correlation function gn{r,N). At the critical point gnir,N) is expected to have the 
following behavior 

—m{N) r 

gnir,N) , (14) 

were the mass gap m(N) vanishes in the infinite volume limit. Surprisingly we 
only see the exponential decay of the spin-spin correlator and not the power fall-off 
underneath it (on a log plot we have a straight line for some range of r). If we 
assume that the inverse mass gap is yet another measure of a characteristic length 
scale for the system, the observed power law dependence is an alternative measure 
of the Hausdorff dimension. 

Looking at the data it is clear that the scaling hypothesis is satisfied as well 
here as for pure gravity. What is surprising is that the extracted values of dn, with 
two exceptions, are almost the same as for pure gravity. The exceptions, for both 
models, are the scaling of the peak height of gun{f^N) and dn obtained from the 
mass gap, both indicating larger values of dn- Why is it that we do not seem to see 
any effects of the back reaction of matter on the fractal dimension? 

A possible explanation would be that the critical region is slightly shifted 
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Ising 

n(r, N) 

dH 5d 


model 

gun [r, N) 

dH 


3-statc Pc 
n{r, N) 
dH x' 


)tts model 

gun{r, N) 
dn X^ 


(a) 

500-1000 

1000-2000 

2000-4000 

4000-8000 

8000-16000 


3.758(53) 2.6 
3.802(55) 0.77 
3.833(56) 1.0 
3.893(61) 0.88 
3.870(87) 0.35 


3.76(12) 0.93 
3.75(15) 1.0 
3.73(12) 2.5 
3.69(09) 3.9 
3.80(10) 0.99 


3.752(63) 0.68 
3.787(65) 0.29 
3.864(63) 1.0 
3.870(73) 0.15 
3.820(97) 0.58 


4.01(26) 2.5 
4.11(18) 1.0 
4.04(22) 3.2 
4.11(19) 0.41 
4.14(15) 0.56 


(b) 

1000-16000 


3.862(74) 1.4 


3.851(53) 4.5 


3.831(32) 2.4 


3.966(64) 12.5 


position 
height 
mass gap 


3.875(53) 
4.01(15) 

4.51 


3.88(19) 
4.36(18) 

(20) 


3.879(29) 
3.900(41) 

4.56 


4.141(58) 
4.424(35)) 

(43) 



Table 3: Extracted values of dn for the Ising and 3-state Potts models coupled to 
gravity. The values are obtained in the same way as for pure gravity (Table 1). 



away from the infinite volume critical coupling at the finite volumes we simulate. 
This is, for example, observed in measurements of the string susceptibility |13|, 
where measured values of 7s peak away from (3c- To check this we have measured 
dfj for the Ising model over an interval of (3. Within errors the extracted value of 
dfj did not change over this interval. 

In the case of pure gravity we see that the scaling of the peak height gives 
better results. If we believe this we get different values for dfj depending on which 
point-point correlator we examine. Looking at n(r, N) we get dn ~ 3.9 for both 
models, and observe no back reaction from the matter. The distribution gun{T,N), 
on the other hand, indicates dn > 4, and indeed gives results that might be consis- 
tent with the values predicted in This is supported by the scaling of the mass 
gap of the spin-spin correlator. We will return to this in the discussion section. 



4 Hausdorff Dimension - Analytic results 

In this section we briefiy review the continuum and matrix model derivations of the 
intrinsic Hausdorff dimension {dn) of the surfaces generated by the coupling of 2d 
gravity to matter ||3|, |^, ^ Q. There are several potentially inequivalent ways to 
define an appropriate measure of the fractal dimensionality of random surfaces. In 
the original paper of [0] two methods were proposed. In the first method one deter- 
mines a power-like relation between two gauge-invariant observables with dimensions 
of volume {V) and length (L) respectively, with dn determined hy V (x L'^" . The 
volume is measured by the cosmological term and the length by the anomalous di- 
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mension of a test fermion which couples to the gravitational field but generates no 
back reaction. This yields 



In the second method one considers the diffusion of a test fermion field and de- 
termines dn by the short-time come-back probability p(r) tx r~^"^'^. The authors 
were able to determine the Hausdorff dimension in a double power series expansion 
in e = D — 2 and — , where D is the classical dimensionality of the surface and c 
is the central charge of the matter coupled to gravity. In |p this second method 
was applied instead to a scalar field - one considers the diffusion equation for a 
random walk on the ensemble of 2d manifolds determined by the Liouville action. 
This yields 

dH = -2^ = 2 V^±V^ , (16) 

where e'^^'^ corresponds to the cosmological constant operator, which has dimension 
one, and e'^-^'i' corresponds to the Liouville dressing of the Laplacian, which requires 
it to be of conformal dimension —1. 

In the matrix-model/dynamical triangulation approach the transfer matrix 
formulation can be used to obtain an expression for the Hausdorff dimension in the 
case of pure gravity ||3|, ^. One finds dn = 4 in agreement with Eq. (^) for c = 0. 

For the case of pure gravity this result can be compared with Using 
matrix model results it is possible to show that 

where p[L\ D) is the number of boundaries separated by geodesic distance D from 
a loop of length L with one marked point, and the scaling variable x = Now 
one can consider the quantities < L" >= dLU^ p{L; D), where a is the lattice 
constant. From Eq. ( |T7|) it can be shown that: 

< L° > ~ const X D^a"^/^ + constDa'^^"^ + constD^ (18) 

< > ~ const X D^a"^/^ + constD"^ (19) 

< L" > ~ const X L>2" (n > 2). (20) 

Then, using the definition (L^) oc r'^H-^ ^ Qj^e can read off the Hausdorff dimension 
dn = 4, which agrees with the continuum result and our numerical results based 
on scaling arguments. This result is not universal because of the explicit lattice 
dependence in < >. One obtains the same result, however, from the second and 
higher moments provided one assumes that < > scales like the area A. The 
result thus appears to be universal. 
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The general situation is, however, far from clear. One case where there is an 
obvious discrepancy seems to be the (2,2fc — 1) series of minimal models coupled 
to gravity. It is possible to extend the continuum Liouville theory analysis to these 
models after taking into account the fact that these non-unitary models possess 
operators in the matter sector with negative conformal dimensions. It is also possible 
to use the results obtained in [15| to calculate the Hausdorff dimensions for models 



(with 'fc' even). We find that the results thus obtained do not agree with each other 
except for the cases k = 1,2. 

The expression for the distribution of loops at a geodesic distance 'D' for the 
(2, 2k — 1) models coupled to gravity (for even 'A;') was computed in |jl^. They find 
that 

p{L;D)^—, 



-x-'^-2(2ct + 1 + x)+ 



ll2r{ar • ^ ' ' T{a + l)_ " ' ^^^^ 

where a = /c — 3/2 and 71, 72 are '/c' dependent constants. Using the same arguments 
as in the case of pure gravity we can compute dn = {2k — l)/{2k — 3) + 1. 

The continuum result of Kawamoto can also be extended to this case, with 
the difference being that the cosmological constant is not the dressing of the identity 
operator but of the operator with the lowest conformal dimension. Similarly the 
dressing condition for the Laplacian is that the Liouville field has dimension — 1 — 

Then one obtains: 

= m ''^ 

a-i-A„. = ^^^^ (2t+l-V32t-15) (23) 

Ak 

dH = , . (24) 

-2k - 1 + ^/32fc^^ 

It is possible to replace the dressing of the Laplacian with the condition that the 
dressing of the Laplacian involves the identity operator and not the minimal dimen- 
sion operator, in which case we obtain: 

~ -2k -1 + ^{Ak^ + 20A: - fj ' ^^^^ 

Thus for this class of models we find an obvious discrepancy between the 
matrix model and the continuum formulations. These models are not, unfortunately, 
amenable to numerical simulations to resolve this disagreement. 

5 Discussion 

We have studied a class of correlation functions defined along geodesic paths in the 
dynamical triangulation formulation of two-dimensional gravity. The critical nature 
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of this theory is revealed in the observation that these correlators satisfy a scaling 
property. The origin of this scaling behavior can be attributed to the existence of 
a dynamically generated length scale in two-dimensional gravity. Furthermore the 
power relation between this linear scale and the total volume allows us to extract 
a fractal dimension characterizing the typical quantum geometry. For pure gravity 
we estimate dn = 3.83(5), which is close to the analytic prediction dn = 4. Our 
numerical method constitutes by far the most reliable method yet investigated for 
extracting this fractal (Hausdorff) dimension. 

Encouraged by this result we have studied two simple spin models coupled 
to quantum gravity — the Ising and 3-state Potts models. As we have indicated 
there are no truly reliable analytic predictions concerning the nature of the fractal 
geometry for these values of the matter central charge. The inclusion of matter fields 
allows us to define two independent correlation functions which we have termed /i 
and f2- The usual geometrical correlator counting the number of sites at geodesic 
distance r is just the sum fi + /2, whilst the weighted difference /i — ^^/2 yields 
the (unnormalized) spin correlator. 

For both types of correlation function in either the Ising or 3-state Potts cases 
we see good evidence for scaling. From the geometrical correlators the Hausdorff 
dimension we extract is statistically consistent with its value for pure gravity. Taken 
at face value this would seem to indicate that the back-reaction of the critical spin 
system on the geometry is insufficiently strong to alter the Hausdorff dimension for 
these values of the central charge. This is supported by our best overall scaling fits 
to the spin correlator, which yield comparable values for dn- 

The picture is somewhat different if we use only the scaling of the peak height 
to estimate a value for dn — now a shift in du is observed to values somewhat above 
four. Indeed these estimates for dn are not inconsistent with the predictions of the 
formula derived in |^]. Since the peak scaling appears to suffer from smaller finite 
size effects than other quantities in the case of pure gravity (it gives dn = 4.040(98)) 
it is possible that it is also a more reliable channel in which to look for signs of back- 
reaction in the case of spin models. These estimates for dn are also favored by 
examining the scaling of the spin correlation length extracted from the normalized 
correlation function. Without good theoretical reasons for believing in such a favored 
channel, however, it is probably more sensible to ascribe the differences in our 
estimates for du to the presence of rather large scaling violations at these lattice 
sizes. 

One alternative scenario might be that the observed effects are due to the 
presence of two linear scales; the geometrical scale and another characterizing the 
critical spin correlations. Thus two fractal dimensions might be possible; one the 
(true) Hausdorff dimension associated with the geometry, and another revealed only 
in the spin channel. We can see how these two scales could coexist by considering the 
correlation functions /i and f2- We have seen that both of these quantities appear to 
scale in a similar fashion. Indeed the exponents associated with the overall volume 
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prefactors are the same. Yet the weighted difference of fi and /2, /i — 
yields the spin correlator which is constrained to have a different dependence on 
the volume. This implies that both fi and ^37/2 are composed of identical leading 
terms together with sub dominant terms. The simplest scenario for fi might be 

/i (r, N) = N-ct> ifi) + Ny^j ifi') {x>y) . (26) 

A similar expression would hold for /2. The idea is that the distribution n{r,N) is 
determined by the lead terms whilst for the precise linear combination making up 

g„„(r, N) this piece cancels and we are left with the subdominant piece -0 ip-')- 

1 

The exponent x is just l — l/dn and the scaling variable fi = r/N''-H = r/lc 
measures the geodesic distance in units of the induced geometrical scale Iq- Similarly 

the exponent y is related to the spin susceptibility exponent y = — and 

1 

/i' = r/N'^s = r /Is is a possible new scaling variable associated with the spin scale 
Is- In flat space the critical spin scale Is is just identified with the linear scale (here 
Ig) and ds = dn- It is not clear on a dynamical lattice that this is necessarily so; 
one could imagine a scenario in which the geometrical scale varies anomalously with 
the spin scale Ig ~ i^s)'^ ■ The quantity uj would then constitute a new exponent 
characterizing the coupled matter-gravity system. 

If this scenario were to be realized then the numerical estimates of these 
exponents would favor a situation in which the spin correlation length diverged more 
slowly with volume than the gravitational (geometrical) scale. This might serve as 
a partial explanation of the observed exponential behavior of the (normalized) spin 
correlator at the critical point - unlike flat space critical models the correlation 
length in a dynamical lattice can never reach the typical linear size of the lattice. 

In the absence of any explicit transfer matrix type solutions for these unitary 
minimal models it would seem that further high resolution numerical work will be 
needed to resolve these important issues. 
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